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ABSTRACT 

We present the discovery of an inhomogenous, low-surface brightness, extended disk-like structure 
around the Andromeda galaxy (M31) based on a large kinematic survey of more than 2800 stars 
with the Keck/DEIMOS multi-object spectrograph. The stellar structure spans radii from 15kpc 
out to ~ 40 kpc, with detections out to R ~ 70 kpc. The constituent stars lag the expected velocity 
of circular orbits in the plane of the M31 disk by ~ 40kms -1 and have a velocity dispersion of 
~ 30kms . The color range on the upper RGB shows a large spread indicative of a population 
with a significant range of metallicity. The mean metallicity of the population, measured from Ca II 
equivalent widths, is [Fe/H] = —0.9 ± 0.2. The morphology of the structure is irregular at large radii, 
and shows a wealth of substructures which must be transitory in nature, and are almost certainly tidal 
debris. The presence of these substructures indicates that the global entity was formed by accretion. 
This extended disk follows smoothly on from the central parts of M31 disk out to ~ 40 kpc with an 
exponential density law of scale-length of 5.1 ± 0.1 kpc, which is similar to that of the bright inner 
disk. However, the population possesses similar kinematic and abundance properties over the entire 
region where it is detected in the survey. We estimate that the structure accounts for approximately 
10% of the total luminosity of the M31 disk, and given the huge scale, contains ~ 30% of the total 
disk angular momentum. This finding indicates that at least some galactic stellar disks are vastly 
larger than previously thought and are formed, at least in their outer regions, primarily by accretion. 
Subject headings: galaxies: individual (M31) — galaxies: structure — galaxies: evolution — Local 
Group 



1. INTRODUCTION 

Within the framework of hierarchical structure forma- 
tion, large spiral galaxies like the Milky Way or An- 
dromeda (M31) arose from the merger of many small 
galaxies and protogalaxies which began coalescing at 
high redshift. During this process, considerable dynami- 
cal and chemical evolution took place, as the primordial 
galaxy fragments were assimilated, stellar populations 
evolved and fed enriched gas back into the interstellar 
medium, and fresh gas was accreted onto the galaxy. The 
disk, the defining component of spiral galaxies, formed 
through dissipative processes, growing from the inside 
out mainly by the accretion of gas with high angular 
momentum with probably only a small component of ac- 
creted stars (Abadi et al. 2003). 

The structure of disks is generally exponential in both 
their radial and vertical directions. A striking feature 
of this component of galaxies is the tendency for it to 
be truncated radially at 3-4 scale-lengths (van der Kruit 
& Searle 1981; Pohlen, Dettmar & Liitticke 2000). This 
distance would correspond to ~ 12 kpc in the Milky Way 
and ~ 20 kpc in M31. Although we still lack a solid 
understanding of how and why these truncations occur, 
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evidence suggests an abrupt change in the star forma- 
tion efficiency across this boundary, perhaps related to 
a shift in the dominant star formation mode (Kennicutt 
1989; Ferguson et al. 1998; Schaye 2004). The fact that 
disk stars are detected beyond the predicted truncation 
radii in the Milky Way and M31, as well as in several 
lower mass galaxies (e.g. NGC 300, Bland-Hawthorn et 
al. 2005; M33, Ferguson et al 2005, in prep), is intrigu- 
ing. Did these stars form in situ? If so, did they form 
a long time ago during the very early stages of disk evo- 
lution, or have they formed much more recently from 
newly-acquired gas? On the other hand, could they have 
been tidally stripped from infalling satellites on close to 
coplanar orbits? 

In order to address these questions and examine galaxy 
formation models in detail, our group has undertaken a 
deep panoramic study of the Andromeda galaxy, which 
when set in contrast with the Milky Way, was to serve 
as a test case against which formation simulations could 
be compared. Studies of the Milky Way are naturally 
blessed with more photons that those of more distant 
systems, and by resolving individual stars and measuring 
their kinematic and chemical properties one can study 
the formation of our Galaxy in a detail and with a mass 
and spatial resolution that can only be dreamed of else- 
where. However, thanks to advances in instrumenta- 
tion, it is now also possible to measure the kinematics 
and chemistry of individual stars in galaxies of the Local 
Group. One of the major benefits that this affords us is a 
global view of a galaxy, freeing us from much of the pro- 
jection uncertainties that our position within the Milky 
Way plagues us with, and yet retaining the fine detail 
that spectroscopic studies of resolved stellar populations 



Fig. 1. — The current coverage of our large panoramic survey of M31 with the INT camera, in standard coordinates (!;,r)). The outer 
ellipse shows a segment of a 55kpc radius ellipse flattened to c/a = 0.6, and the major and minor axis are indicated with straight lines out 
to this ellipse. A DSS image of the central galaxy has been inserted to scale. The map shows the distribution of RGB stars throughout the 
halo. The most striking structure in the map is a vast flattened inner structure, of major axis extent ~ 4 degrees, a messy inhomogenous 
entity which envelops the familiar bulge and disk of M31. 



provide. 

Located at a distance of 785 ± 25kpc ((m — M)q = 
24.47 ± 0.07; McConnachie et al. 2005), Andromeda is 
the closest giant spiral galaxy to our own, and the only 
other giant galaxy in the Local Group. As such it is 
often considered to be the "sister" galaxy to the Milky 
Way. Recent studies suggest that the two galaxies do in- 
deed have very similar total masses (including the dark 
matter, Evans et al. 2000; Ibata et al. 2004), yet there 
are significant differences between them. M31 is slightly 
more luminous than the Milky Way, it has a higher ro- 
tation curve, and a bulge with higher velocity disper- 
sion. M31 possesses and a globular cluster system with 
~ 500 members, approximately three times more nu- 



merous than that of the Milky Way, and contains sev- 
eral members of a type of luminous extended cluster not 
present around the Galaxy (Huxor et al. 2005). The 
disk of Andromeda is also much more extensive, with 
a scale-length of 5.9 ± 0.3 kpc (R-band value corrected 
for a distance of 780 kpc, Walterbos & Kennicutt 1988) 
compared to 2.3 ±0.1 for the Milky Way (Ruphy et al. 
1996); but which is currently forming stars at a lower 
rate than the Galaxy (Avila-Reese, Firmani & Hernandez 
2002; Walterbos & Braun 1994). Though possibly the 
consequence of low-number statistics, it is tempting to 
attribute significance to the fact that Andromeda has 
four dwarf elliptical galaxies (M32, NGC 205, NGC 147, 
NGC 185) among its entourage of satellites and no star- 
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forming dwarf irregulars, whereas the Milky Way has 
none of the former but two of the latter. However, it 
is perhaps in their halo populations that the differences 
between the two galaxies are most curious and most in- 
teresting. It has been known for many years that M31 
does not possess the "text-book" stellar halo model de- 
veloped from observations of the Milky Way. Minor 
axis halo fields have shown startlingly high metallicities 
([Fe/H] ~ —0.5) and an i? 1 / 4 -law density profile (Dur- 
rell, Harris & Pritchet 2004), unlike the R~ 35 behavior 
of the metal-poor Galactic halo (Chiba & Beers 2000). 
Nevertheless, recent studies of M31 confirm that an old, 
metal-poor halo lurks beneath (Holland 1996; Durrell, 
Harris & Pritchet 2001; Reitzel & Guhathakurta 2002; 
Brown et al. 2003). Both the de Vaucoulcurs profile and 
the high mctallicity are suggestive of an active merger 
history at the time of halo (or bulge) formation. 

The photometric dataset that forms the foundation of 
this study was obtained with the Wide Field Camera at 
the 2.5m Isaac Newton Telescope; the survey and pre- 
liminary results were presented in Ibata et al. (2001b) 
and Ferguson et al. (2002), though the full completed 
survey will be discussed in detail in a forthcoming ar- 
ticle (Irwin et al., in prep). The density of Red Giant 
Branch (RGB) stars throughout the survey region, which 
now covers ~ 40 sq. degrees, is reproduced in Figure 1, 
onto which we have superimposed the familiar Digital 
Sky Survey image. The solid angle covered by the sur- 
vey corresponds to a projected area of ~ 7400 kpc 2 at 
the distance of M31. As shown in the earlier publica- 
tions in this series, the INT survey revealed a wealth of 
substructure throughout the "halo" of M31. However, 
the most prominent of the details that we uncovered was 
the presence of a gigantic extended flattened structure 
around M31, spanning up to ~ 40 kpc from the galaxy 
centre, almost twice the maximum extent of the bright 
inner disk. Due to the expectations from galaxy forma- 
tion simulations with cold dark matter, our initial inter- 
pretation of this structure favored the possibility that it 
was a flattened inner halo (Katz & Gunn 1991; Summers 
1993; Steinmetz & Miiller 1995) with significant spatial 
sub-structure (Klypin et al. 1999; Moore et al. 1999). 
However, the proximity of M31 has allowed us to follow 
up these detections with deep Keck/DEIMOS spectro- 
scopic observations (initial results of the spectroscopic 
survey were presented in Ibata et al. 2004). It is the ra- 
dial velocities and metallicities derived from these deep 
spectroscopic data that we present and analyze in this 
contribution. 

The layout of this paper is as follows. In §2 we first 
discuss the global geometry of M31, the inclination angle 
of the disk and fit an exponential function to the major 
axis photometry. The spectroscopic observations are pre- 
sented in §3, and the resulting kinematics in §4. In §5 
we examine several fields in more detail, presenting also 
their mean metallicity and fitting global trends. These 
data are interpreted in §6, where we present several pos- 
sibilities for the formation of the structure, and we finish 
in §7 with the conclusions of the study. 

2. SIMPLE GLOBAL STRUCTURAL MODEL OF M31 

Before describing the observations and analysis of the 
kinematic dataset, it is useful to discuss the large-scale 
structure of M31. Andromeda is observed highly inclined 



to the line of sight with i = 11° (Walterbos & Kennicutt 
1988), and it is seen on the sky with a position angle of 
38° East of North (Walterbos & Kennicutt 1987). Within 
an angular distance of 90' (20.4 kpc), Walterbos & Ken- 
nicutt (1988) find a good fit with the sum of a bulge 
model together with a disk model of R-band scale-length 
5.9 ± 0.3 kpc (corrected for a distance of 780 kpc). 

To compare these findings to our INT survey, we mea- 
sure the source density in a color-magnitude box designed 
to select the brightest RGB stars 1.5 magnitudes below 
the RGB tip. Sources with I-band magnitudes between 
20.0 < io < 21.5 and colors 1.0 < (V - i) < 4.0 are se- 
lected. (The Schlegel et al. 1998 dust maps are used 
throughout this work to correct the INT photometry for 
extinction). We correct the star-counts for crowding us- 
ing the approach of Irwin & Trimble (1984), and subtract 
a small background count as determined from fields at 
the edge of the survey. 

The de-projected positions of the RGB stars are binned 
in a polar coordinate grid with 5° azimuthal bins and 
0.5 kpc radial bins. The azimuthal counts at a given 
radius are then medianed together in groups of 5 (i.e. 
covering 25°), and we step around the circle in intervals 
of 15°. The upper panel of Figure 2 shows the 4 pro- 
files within a de-projected angle of ±35° of the major 
axis (i.e. ±7.9° projected), assuming an inclination of 
i = 11°, with the additional thick-lined profile showing 
the median of all 14 bins of 5° width. Positive radius 
means that the medianed area is on the North-eastern 
major axis, and negative radius on the South-western 
major axis. The abscissa R is a de-projected radius, that 
assumes the stated inclination. The artificial drop in the 
source counts inside of R ~ 18 kpc is due to crowding. 

Evidently the major axis profile is close to exponential; 
it is fit in the lower panel with a model (dashed line) of 
scale- length 5.1 ± 0.1 kpc. This scale- length is similar to 
that of the inner disk, so it appears that the large-scale 
properties of the outer disk structure follow on smoothly 
from those of the inner disk. 

The upper panel of Figure 3 shows the profiles for all 
azimuthal bins. The large discrepancies between the dif- 
ferent profiles arc very obvious, with the less steep pro- 
files corresponding to azimuthal bins towards the minor 
axis. The higher number densities on the minor axis are 
largely due to the huge bulge component in M31. In- 
deed, in the bulge-disk decomposition of Walterbos & 
Kennicutt (1988), the bulge dominates over the disk at 
all radii on the minor axis. However, part of this varia- 
tion may also come from a change in the inclination of 
the galaxy at large radius. To explore this possibility, 
we fitted the value of i by assuming that the most likely 
inclination is that that has the lowest root mean square 
dispersion about the median profile in the distance in- 
terval 20 < R < 40 kpc. The corresponding profiles are 
shown on the lower panel of Figure 3, together with the 
medianed profile (thick black line) . The fitted inclination 
angle is i = 64.7°. 

With this smaller inclination, the surface density pro- 
files are now much more similar, yet there are still vari- 
ations between the profiles of up to a factor of ~ 2- 
3 at a given radius, due in part to the presence of 
obvious "lumps". Fitting the median profile between 
20 < R' < 40 kpc with an exponential-law yields a 
scale-length of 7.8 ± 0.07 kpc. This larger scale-length 
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Fig. 2. — The top panel shows the density profiles close to the 
major axis. Positive abscissa values correspond to de-projected 
radii towards the North-eastern side of M31. The thin lines are 
the density profiles medianed over 5 5°-de-projected azimuthal 
bins, assuming a galaxy inclination of i = 77° . The 4 profiles 
to the North-east cover the de-projected angles —35° < 9 < —10°, 
-20° < 9 < 5°, -5° < 9 < 20°, 10° < 9 < 35°, with the 4 pro- 
files to the South-west being identical cuts rotated by 180°. The 
thick black line corresponds to the median of all 14 azimuthal bins 
between —35° < 9 < 35°. The lower panel shows an exponen- 
tial fit to the medianed profile between 20 < R < 40 kpc, with a 
scale-length 5.1 ± 0.1 kpc. 
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Fig. 3. — The top panel shows medianed density profiles at 
different position angles around M31, offset at 15° from each other, 
assuming a galaxy inclination of i = 77°. As in Figure 2, each 
profile is constructed from a median of 5 consecutive azimuthal 
bins of width 5° . The thick black line corresponds to the median of 
all azimuthal bins. The lower panel shows similar profiles, but for a 
an inclination angle of i = 64.7°, fitted to minimize the dispersion 
around the median profile of all azimuthal bins over the interval 
20 < R' < 40 kpc. 



value compared to the previous major axis exponential fit 
demonstrates a failure in our simple tilted disk model. A 
more complex model with a bulge component is required 
to fit the central regions of the galaxy more reliably, yet 
this cannot be constructed from our INT data which are 
highly incomplete inside of R ~ 15 kpc. 

In the following discussion, we will present results 
based on inclination values of % = 77° and i — 64.7°; 
this will help clarify the uncertainty due to a possible 
tilt of the outer disk with respect to the inner disk. De- 
projected quantities, such as the de-projected radius R 
or the velocity lag vi ag (defined below) will implicitly as- 
sume an inclination of i = 77°; primed versions of these 
variables (e.g. R' , v' lag ) will instead correspond to the 
lower inclination case of i = 64.7°. 

3. OBSERVATIONS 

During eight Keck nights (~ 75% useful time) in Sept. 
2002, 2003, and 2004, we undertook a DEIMOS program 
to spectroscopically follow-up the detected substructures 
in M31 to measure their kinematics, and to begin a com- 
plete kinematic portrait of M31 targeting systematically 



disk and halo fields. The observing runs have employed 
the standard high- res DEEP2 slitmask approach (e.g., 
Davis et al. 2003) covering the wavelength range 6400- 

o o 

9000 A, with a spectral resolution of ~0.6 A. The ob- 
servations have also pioneered two new approaches with 

DEIMOS: a band-limiting ~ 300 A Ca II triplet filter to 
multiplex ~ 4" slitlets in the spectral direction yielding 
as many as 800 slits per mask, and a 'fibre-hole' 0"7 slit- 
let approach packing ~ 600 holes per mask. This 'fibre- 
hole' approach adopted in Sept. 2004 proved to be very 
successful, giving Poisson-limited sky-subtraction (down 
to i=21.5) by assigning holes to monitor the sky spec- 
trum; longer slitlets simply waste valuable CCD real- 
estate in regions of high target density. The fiber-hole 
diameters were milled at 0"7 to match the median see- 
ing. The dataset presented in this paper is the sum of 
all the 'standard mask' spectra plus 'fibre- holes' spectra 
obtained up to September 2004 that do not use the band- 
limiting filter. A total of 28 fields were observed in this 
way. The layout of these spectroscopic fields is displayed 
in Figure 4: the 5 fields D1-D5 probe the M31 disk, 
the 16 fields F1-F16 probe both the halo and what we 
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Fig. 4. — The distribution of Keck/DEIMOS fields in relation to 
the survey region depicted in Figure 1 . For clarity, we have divided 
these fields into three groups. The fields in the giant stream are 
prefixed by the letter S, and we keep the same field names as in 
Ibata ct al. (2004), with the addition of fields S24, S26, and S27. 
Those fields located within the inner ellipse arc denoted Dl - D5. 
The remaining fields are named Fl - F16. The small square labeled 
"Brown" on the South-eastern minor axis marks the location of the 
deep photometric study of Brown et al. (2003), whose relevance to 
the present work we discuss in §6. 



considered to be "halo" overdensities, while 7 fields were 
assigned to investigate the dynamics of the giant stream 
(Ibata et al. 2001b) which is seen clearly in Figure 1 as 
an elongated structure to the South-east of M31. 

For each of these fields except the stream fields, the 
highest priority targets were selected for observation 
from the INT survey within a color-magnitude box with 
I-band magnitudes between 20.5 < i < 22.0 (without red- 
dening correction) and colors 1.0 < (V — i)o < 4.0. This 
broad selection was designed to pick out a wide range of 
red-giant branch stars, containing both metal-poor and 
metal-rich populations (further details will be provided in 
Irwin et al. 2005, in prep). Other non-saturated targets 
brighter than i = 22 were also chosen by an automated 
selection algorithm at lower priority to fill in available 
space on the spectrograph detector. 

Each mask was observed for typically three integra- 
tions of 20min each. The spectroscopic images were 
processed and combined using the pipeline software de- 



veloped by our group. This software debiasses, per- 
forms a flat-field, extracts, wavelength-calibrates and 
sky-subtracts the spectra. The radial velocities were 
measured by cross-correlating against a model template 
with Gaussian functions at the rest-frame wavelength po- 
sitions of the Ca II triplet lines (the technique used is sim- 
ilar to that discussed by Wilkinson et al. 2004). An esti- 
mate of the radial velocity uncertainty is then obtained 
for each star from the dispersion of the radial velocity 
measured from the three Ca II triplet lines separately. 
The accuracy of these data are very good for such faint 
stars, with typical uncertainties of 5-10 km s -1 . The 
sample presented in this paper consists of 2834 stars with 
velocity uncertainties < 20 km s _1 that fall into the above 
RGB color-magnitude selection window. 

4. RESULTS 

The velocity distribution of all the targeted stars with 
(V — i)o > 1-0 is displayed on the top panel of Figure 5. 
This distribution is a mixture of bona-fide M31 stars 
plus a contamination from Galactic stars. We have su- 
perimposed the expected Galactic contribution over a 
0.25° area in the direction of M31 according to the Be- 
sancon Galactic starcounts model (Robin et al. 2004), 
onto which we have applied the color-magnitude selection 
window used for spectroscopic target selection. Since we 
have not observed every star in each field, and since the 
ratio of M31 to Milky Way stars varies dramatically from 
the inner fields to the outer fields, we cannot easily scale 
the model predictions. For the illustrative purposes of 
Figure 5, we have chosen a normalization to match the 
Besancon model to our observed distribution for veloci- 
ties higher than the peak in the model. We will return 
to the issue of Galactic contamination in our analysis of 
field F16 (Figure 24), placing an upper limit on the model 
normalisation, which shows that Galactic contamination 
is low at Heliocentric radial velocities below —100 kins" 1 . 
It is clear from a visual comparison of the distributions 
in Figure 5 that we observe a large number of stars in 
significant excess over the expected distribution of con- 
taminants. The apparent mismatch between the peak 
in the observed distribution and the Besangon model is 
due to the large number of stars observed in the disk 
fields D1-D5 on the North-eastern side of Andromeda, 
which combine in this summed-up plot to give a broader 
peak centered on — 60kms _1 . If we instead sum up the 
11 fields on the South-western side of M31 (fields Fl- 
F4 and all the stream fields), the rotation of the M31 
disk separates Galactic stars from the majority of stars 
in Andromeda, allowing a cleaner discrimination. This 
is shown on the bottom panel of Figure 5. These con- 
siderations show that by selecting stars with velocities 
< — 100 km s -1 we can drastically minimize the Galactic 
contamination, incurring only the penalty of including a 
small fraction of Galactic halo stars into our sample; fur- 
thermore these Galactic halo stars are expected to have 
a relatively flat distribution in radial velocity. 

Due to the wide coverage of the spectroscopic fields, 
it is useful to compare the observed velocity of a star 
with the expected radial velocity of a population moving 
in circular orbits. We construct a simple model assum- 
ing the geometrical model discussed in §2. We adopt 
the circular velocity curve compilation of Klypin, Zhao 
& Somerville (2002), which joins CO measurements from 
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Fig. 5. — The top panel shows the velocity distribution of the 
2834 RGB star candidates in all fields shown in Figure 4 that pass 
our color-magnitude and spectrum quality criteria. The green-lined 
histogram shows the prediction of the Besangon Galactic model 
for an area of 240 sq. arcmin within our color-magnitude selection 
window. With this normalization, the model fits the high-velocity 
edge of the observed velocity distribution. The bottom panel shows 
the summed distribution of stars in the fields F1-F4 plus all the 
stream fields. 



Loinard, Allen & Lequeux (1995) for R < 10 kpc with H I 
measurements by Brinks & Burton (1984) at larger radii; 
these data are reproduced in Figure 6, along with our 
spline fit which is shown adjusted for the small inclina- 
tion correction. The dip in the rotation curve near 5 kpc 
is of course an artifact of non-circular motions in the 
gas at that location, and is unimportant for the present 
work since our fields all lie beyond that radius. The cor- 
responding map of the projected rotation velocity over 
the inner region of our survey is shown in Figure 7. 

In reality, the M31 disk has finite thickness, and due to 
the highly inclined orientation that it presents to us, lines 
of sight through the disk will probe a substantial depth. 
This is a concern as the simple circular velocity model of 
Figure 7 could be misleading, since stars at a given (£, rj) 
on the sky can be located at different radii in M31, and 
will also present different projections to us. In Figure 8 
we show the velocity uncertainty caused by a disk scale- 
height of 0.35 kpc, a value comparable to the scale-height 
of the old disk of the Milky Way (Gilmore & Reid 1983). 
For most of the survey area the additional uncertainty to 
the circular velocity model due to a disk of that thickness 
is less than lOkms^ 1 , though certain regions within the 
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Fig. 6. — The rotation curve of M31. The points show a rotation 
curve compilation by Klypin, Zhao & Somerville (2002). which is 
derived from CO rotation velocities measurements (Loinard, Allen 
& Lequeux 1995) for R < 10 kpc, and from the rotation velocity 
of the H I (Brinks & Burton 1984) beyond that radius. The dot- 
ted line shows a spline fit to these data, while the continuous line 
indicates the correction for the 77° inclination of the inner disk of 
this galaxy. 
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Fig. 7. — A simple model for the expected velocity of stars 
on circular orbits as a function of position around M31, made by 
projecting the spline fit of Figure 6 using the geometrical model 
proposed in §2. For comparison with Figures 1 and 4, we have 
superimposed the inner ellipse, and a section of the outer ellipse 
displayed in those figures onto this diagram. The radial velocity 
contours are shown at intervals of 50kms — 1 . 



inner ellipse the uncertainty is > 30kms . We will bear 
this consideration in mind when discussing the kinematic 
properties of the disk. 



4.1. Inner fields 
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Fig. 8. — Since the Andromeda galaxy is observed close to edge- 
on, lines of sight through the disk can transverse a significant depth, 
unless the disk is infinitesimally thin. Here we illustrate the uncer- 
tainty caused by a 350 pc disk scale-height on the circular velocity 
predictions displayed in Figure 7. The contours show the levels 
where the uncertainty in the circular velocity model due to this 
disk thickness attains values of (from the outside in) of lOkms , 
20kms~ 1 and 30kms — 1 . We include this additional error in when 
calculating the lag in velocity of stars behind the circular speed. 
However, over most of the survey region this additional uncertainty 
is less than lOkms -1 . 



We now proceed to present the velocity data. The 
upper panel of Figure 9 shows the (Heliocentric) ra- 
dial velocity measurements as a function of a;, the ma- 
jor axis distance 7 for the disk fields labeled D1-D5 and 
which lie within the inner ellipse of Figures 1 and 4 
(corresponding to a disk radius R < 27kpc); a small 
number of stars from field F13 also fall into this selec- 
tion. The dashed line here marks a Heliocentric veloc- 
ity of v = — 300kms _1 , the adopted systemic velocity 
of M31. Non-rotating populations such as might be ex- 
pected for the stellar halo, will scatter in a broad distri- 
bution around this line. In these North-eastern fields, the 
disk lies at radial velocities above the systemic value, so 
counter-rotating populations can be identified from their 
location below the dashed line. The bottom panel shows 
the result of subtracting the model of Figure 7 from the 
velocity data on a star by star basis (rather than using 
a single value for a whole field). Stars on circular orbits 
in the plane of the disk lie at a disk-corrected velocity of 
v iag = (emphasized with a dashed line). Field Dl (at a 
major axis distance of x ~ 7kpc) appears to be peculiar 
with an intriguing kinematic sub-structure at Heliocen- 
tric velocity of v ~ — 250 km s -1 ; we will discuss this 
structure in a subsequent contribution. It is interesting 
to note that this sub-structure was only detected through 

7 We define the major axis distance as the distance from the 
minor axis towards the North-east parallel to the major axis. 
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Fig. 9. — Radial velocity as a function of major axis distance for 
fields D1-D5 (with a few stars from F13 that lie within the inner 
ellipse of Figure 4). The upper panel shows the Heliocentric radial 
velocity of the stars, with the dashed line marking the systemic ve- 
locity of M31. The bottom panel shows the same data "corrected" 
for the expected disk rotation given by the model shown in Fig- 
ure 7; now the dashed line marks the location of stars on circular 
orbits. Note the obvious structure that is present in the field at 
major axis distance of ~ 7kpc (field Dl). 



its kinematic signature, as there is no clear overdensity 
in our maps of RGB stars in this field. 

However, the disk fields D2-D5 (located at major axis 
distances between 12 kpc < x < 23 kpc) show a very 
similar distribution of stars, with a broad halo compo- 
nent below vi a g = — 150kms , and a well-populated 
peak lagging circular orbits only slightly. In Figure 10 
we show the distribution of disk lag velocities for fields 
D2-D5, together with a Gaussian fit of the data between 
— 150 kms^ 1 < vi a g < 100 kms -1 , which has a mean of 
viag = — 34kms~ 1 with dispersion of a v = 47kms _1 
i v 'iag = — 18kms _1 , av = 47kms~ 1 ). The dash-dotted 
line shows a similar maximum-likelihood Gaussian fit, 
this time taking into account the estimated velocity un- 
certainty on the individual velocity measurements, and 
also the uncertainty in the disk rotation model (from 
Figure 8). This fit, which has vi ag = — 37kms _1 and 
a v = 44kms~ 1 (v'lag = — 23kms~ 1 , ov = 45kms~ 1 ), 
shows that the dispersion of the population is not arti- 
ficially enhanced to any significant extent by our line of 
sight depth through these fields. 

As discussed above, contamination from Milky Way 
stars can be problematic for v > — lOOkms ; however 
the M31 disk has very high contrast in these inner fields 
(D1-D5), with a source surface density in our CMD se- 
lection window that is 20 times higher than in the halo 
field F7, so we expect the Galactic contamination to be 
small, 5% at the very most. Thus we find that the dom- 
inant population in the inner regions of M31 is a thin 
disk population with a dispersion of a v = 44kms _1 
that trails circular orbits by Vi ag = — 37kms~ 1 ; these 
values are reasonably similar to the measurements of 
the velocity dispersion of the Galactic disk in the So- 
lar Neighborhood where the velocity ellipsoid has axes 
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Fig. 10. — The summed distribution of velocity lags of fields 
D2-D5 (i.e. stars with major axis distances between 12kpc to 
23kpc in Figure 9). The large peak between — 150kms _1 < 
v lag < 100 km s -1 can be fit with a Gaussian (dashed line) of 
mean vi ag = — 34kms~ 1 and dispersion <j„ = 47 kms" 1 (v'[ ag 

■ a v i = 47kms~ 1 ). The fit does not change signifi- 

cantly if we take into account the estimated velocity errors and 
the disk rotation model uncertainties displayed in Figure 8, we ob- 
tain then (dot-dashed line) vi ag = — 37kms" 1 , a v = 44kms — 1 
( v 'lag = —23 kms" 1 . a v i = 45kms~ 1 ). 
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Fig. 11.— As Figure 9, but for the 16 fields F1-F16 that lie 
outside of the inner ellipse shown in Figures 1 and 4, and which 
no not lie in the giant stream. In most of these fields the domi- 
nant population of stars is seen with velocities lagging the circular 
velocity by velocities between Okms -1 to ~ 100 kms -1 . 



(<Jr, <J<f,, a z ) = (38 : 24 : 17) kms 1 and the population 
lags the circular velocity with an asymmetric drift of 
~ 20 kms" 1 (Dehnen & Binney 1998). The asymmet- 
ric drift of the M31 population will be discussed further 
in §6.1. Note that the very inclined orientation of M31 
means that we do not have access to the vertical disper- 
sion o z . The measured a v is a mixture of the radial and 
azimuthal components ctr and o^, dependent on the field 
under study. For simplicity we neglect this distinction in 
the present discussion, referring to a single a v . 



Fig. 12. — As Figure 11, but showing the velocity as a function 
of minor axis distance. 



4.2. Outer fields 

Having shown that the inner fields have kinematics 
typical of a canonical thin disk, we now examine the 
outer field data. The fields lying on the giant stellar 
stream (those prefixed by the letter 'S' in Figure 4) are 
removed from the sample for the sake of simplicity. This 
leaves a total of 16 fields marked F1-F16 on Figure 4. 

The measured radial velocities in these fields are dis- 
played as a function of major axis distance on the top 
panel of Figure 11. For clarity, we have removed from 
this sample all stars with heliocentric velocities v > 
— 100 kms -1 , thereby minimizing the Galactic contam- 
ination (c/. Galactic model in Figure 5). There is a 
wealth of kinematic substructure in these fields (Chap- 
man et al. 2005, in prep), and a significant halo popula- 
tion seen scattered around the systemic velocity (dashed 
line); however the dominant population in these fields 
again follows closely the disk rotation. This can be seen 
on the bottom panel of Figure 11, where we have cor- 
rected the data on the upper panel for the rotation model 
of Figure 7. In the fields at major axis distance be- 
tween —20 kpc < x < —30 kpc, the numerous popula- 
tion stands out clearly with "disk-corrected velocity" be- 
tween Okms -1 and 100 kms -1 . On this South-western 
side of M31, these velocities correspond to a lagging be- 
hind the H I. On the opposite side of the disk, at major 
axis distances x > 5 kpc, a similar concentration is found 
with "disk-corrected velocity" between —100 kms -1 and 
Okms -1 , this again corresponds to a lagging velocity be- 
hind the H I gas. 

For completeness, in Figure 12 we show the same ve- 
locity data as a function of minor axis distance 8 , y. It 
is clear from this diagram that the halo fields F5 and 
F7 (the two fields at minor axis distance y < —15 kpc) 
do not contribute significantly (if at all) to the peaked 
velocity structure that rotates as a disk. 

Figure 13 shows the distribution of velocity lag vi ag in 
the 16 fields under consideration. The velocity lag is de- 

8 We define the minor axis distance as the distance from the 
major axis towards the North-west parallel to the minor axis. 
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Fig. 13. — Distributions of the velocity lag behind the circular 
velocity model for the 16 fields displayed in Figure 11. In each panel 
the three numbers show, respectively, from top to bottom, the 
major axis distance, the minor axis distance, and the de-projcctcd 
radius R of the fields. The field name is also indicated. 



fined to be equal to the disk-corrected velocity for fields 
at positive major axis distance, but is inverted for fields 
at negative major axis distance. Each panel is labeled 
with the field name and as well as three numbers which 
give the major axis distance x, the minor axis distance y 
and the de-projected M31-plane radius R, of the centre 
of each field, all in kpc. A striking feature of these distri- 
butions is the presence of a narrow component in most 
of the fields, which lags the circular velocity by merely 
~ 50 kms -1 , that is to say a population which has dis- 
tinctly disk-like kinematics. Yet this is a startling find- 
ing, since these fields are extremely distant from M31, 
lying at radial distances between 30 kpc and more than 
50 kpc. That this behavior is the most salient kinematic 
feature of the survey is demonstrated in Figure 14, which 
shows the distribution of velocity lag as a function of the 
M31-plane radial distance, R, for all stars in the survey 
in non-stream fields with R < 80 kpc. The majority of 
the stars are found with disk lag velocity in the range 
— 100 kms" 1 < vi a g < Okms -1 , and therefore take part 
in a large-scale rotational motion. 

Summing the velocity lag distributions in all the outer 
non-stream fields (except for F6 to remove any concern 
regarding contamination from NGC 205 — this will be 
discussed further in §5.4) gives the distribution shown in 
Figure 15. This distribution is very similar to that pre- 
viously constructed for the inner fields (cf. Figure 10), 
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Fig. 14. — The distribution of velocity lags as a function of M31- 
plane de-projected radial distance in all non-stream survey fields 
at R < 80 kpc. The majority of these stars display small velocity 
lags (—100 kms -1 < Vi ag < Okms -1 ), and are therefore rotating 
rapidly. The dashed line shows the location of circularly-rotating 
stars. 
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Fig. 15. — The summed distribution of velocity lags in all 
fields in Figure 13, except F6. A Gaussian fit to the data be- 
tween — 150kms — 1 < vi ag < lOOkrns - 1 is superimposed (dashed 
line), which has mean at vi ag = — 54kms — 1 and a dispersion of 
a v = 37kms — 1 (v'lug = — 41 kms -1 , <j v i = 38 kms -1 ). 



showing again a broad halo distribution with the ad- 
dition of a very strong peak, in this case centered at 
viag = — 54kms" 1 with a dispersion of a v = 37kms -1 
(v'lag = —41 kms" 1 , a v > — 38 kms" 1 ). 

These data unambiguously show that the bulk of the 
stars in our survey fields share a common kinematic prop- 
erty: they co-rotate with the H I gas, with velocities that 
lag that of the gas by ~ 40 km s -1 , and possess a velocity 
dispersion of ~ 40 kms" 1 . There therefore appears to be 
a vast stellar disk-like structure around the Andromeda 
galaxy, distributed on an unprecedented scale. 

5. OVERVIEW OF SELECTED OUTER FIELDS 

We next present in detail the data from 7 of the outer 
fields, selected from Figure 13 for having a particularly 
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strong disk-like peak in their velocity distributions. The 
selected fields are Fl, F2, F3, F6, F8, F12 and F16. 

5.1. Field Fl 

Field Fl lies at the South-western boundary of the 
conspicuous dark and messy structure in Figure 1 that 
surrounds the bright inner disk seen in photographic 
surveys. The most massive of the M31 outer globular- 
clusters, named "Gl" also lies in this field, and we re- 
fer to the structure at this location as the Gl clump. 
We summarize the structural, kinematic and photomet- 
ric data on this field in Figure 16. The upper left- 
hand panel gives the location of the field with respect 
to the large survey area, which wc zoom into in the 
upper right hand panel. The bottom left panel shows 
the velocity distribution of the stars in our RGB se- 
lection window, together with a Gaussian fit centered 
at viag = — 69kms~ 1 with dispersion a v = 30kms -1 
{ v 'iag = — 51kms _1 , cr v i = 29kms~ 1 ). The bottom 
right-hand panel shows the CMD derived from the INT 
survey over the 16.7' x 5' field of view of the DEIMOS 
instrument. Those stars with velocities within 2a of the 
mean of the Gaussian fit are marked with a filled circle. 
The four curves superimposed on this CMD are globu- 
lar cluster fiducials, with metallicities, from left to right 
of [Fe/H] = -1.91, [Fe/H] = -1.29, [Fe/H] = -0.71 and 
[Fe/H] = —0.2. The stars that partake in the strong 
peak evidently scatter in a wide distribution around 
[Fe/H] ~ —0.7. Both the kinematics and the spectro- 
scopic metallicity (presented in §5.9) that we derive in 
this field are in excellent agreement with the findings of 
Reitzel, Guhathakurta & Rich (2004). 

The spatial structure of this grouping of stars is very 
interesting. An inspection of the top right-hand panel 
(and also Figure 1) shows that this structure is an irreg- 
ular lump - somewhat elongated E-W and apparently 
distinct from the disk. The structure is very large, be- 
ing more than 10 kpc long. It is highly unlikely that this 
structure is bound, as this would require some ~ 10 9 M© 
of dark matter, given the velocity dispersion and extent. 
Instead, the disorderly morphology argues in favor of this 
being tidal debris. 

5.2. Field F2 

Field F2 (Figure 17) is located close to Field Fl, but 
in a prominence that descends southwards of the disk. 
The velocity distribution in field F2 is similar to that of 
Field Fl, but contains an additional narrow component, 
which is seen as a spike near vi ag = in the bottom left 
hand panel of Figure 17. Fitting two Gaussian functions 
yields a narrow peak at viag — — 4kms~ 1 with disper- 
sion a v — 7kms~ 1 {v'i ag = llkms -1 , a v < — 5kms _1 ), 
and a broader peak with viag = — 46kms~ 1 and a v = 
SSkms" 1 (Viag = -29kms- 1 , a v > = SSkms" 1 ). The 
stars that gives rise to the narrow kinematic spike are 
dispersed spatially over the whole DEIMOS field, so the 
structure is not simply a cluster. The small velocity dis- 
persion and the mean velocity almost exactly on the ex- 
pected circular rotation imply that the the population 
is almost certainly young. The population in the broad 
kinematic peak is again very similar to the peak that we 
observe in field Fl and in the inner fields D2-D5. 
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Fig. 16. — The upper left panel shows the location of the field 
Fl in our survey region (compare to Figure 1), with a 1° X 1° 
box marking the location of the zoomed-in view displayed on the 
right-hand panel, which is centered on the centre of this DEIMOS 
field. The white holes present in this map are due to the halos 
of bright stars. The bottom left-hand panel shows the distribu- 
tion of the velocity lag behind the circular velocity in this field, 
together with a single-component Gaussian fit, which has a mean 
of viag = — 69kms _1 and dispersion of a v = 30kms _1 (v'i ag = 
— 51kms — , <j v i = 29kms — 1 ). The bottom-right panel shows the 
color-magnitude diagram of the field, with small dots showing the 
full INT photometry within the area of the DEIMOS field of view. 
The filled circles are the stars with observed velocities within ±2<r 
of the fitted Gaussian on the bottom left-panel. For comparison, 
we have superimposed fiducial cluster sequences for (from left to 
right) NGC 6397, NGC 1851, 47 Tuc and NGC 6553, which have 
metallicities of [Fe/H] = -1.91, [Fe/H] = -1.29, [Fe/H] = -0.71 
and [Fe/H] = —0.2, respectively. 



5.3. Field F3 

Field F3 is located towards the South-west of M31, on 
a dense elongated structure that trails behind the disk 
in a way that is reminiscent of a tidal stream or the 
end of a spiral arm (see Figure 18). Though the velocity 
dispersion of a v = 15 kms -1 is relatively low in this field, 
the color-magnitude position of the stars that populate 
the kinematic peak is very similar to that seen in other 
fields. 

5.4. Field F6 

Figure 19 shows the velocity and color-magnitude 
data for field F6, centered on the Northern edge of 
NGC 205, one of the dwarf elliptical satellite galaxies 
of Andromeda. The distribution of velocities has a nar- 
row peak centered at vi ag = 2kms^ 1 with dispersion 
<j v = 36kms _1 (v'iag — 4kms , a v < = 37kms _1 ). Sur- 
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Fig. 17. — As Figure 16, but for field F2. In this case we 
have fit a two-component Gaussian model: we find a narrow peak 
centered at vi ag = — 4kms" 1 with dispersion of a v = 7kms — 1 

( v ' lag = llkms - 1 , a v i = 5kms _1 ), and a broader peak with 
viag = — 46kms — 1 and a v = 35kms _1 (y' i ag = — 29kms — 1 , 
<j v i = 33 kms -1 ). 



Fig. 18. — As Figure 16, but for field F3. The Gaussian model 
has vi ag = — 53kms _1 and a v = 15kms — 1 (v'i ag = — 37kms — 1 , 
a v i = 15kms -1 ). 



cr v i = 35 kms 1 ). 



prisingly, this distribution is very similar to what has 
been found previously in other fields: the velocity dis- 
persion is small and consistent with a circularly rotating 
population. Figure 20 displays the velocity and spatial 
data in more detail, with the middle panel showing the 
target positions with respect to a sketch of NGC 205. 
The upper panel gives the heliocentric velocity measure- 
ments, while the bottom panel shows the velocity lag as 
a function of standard coordinate £. It is interesting to 
note that if some stars in this field were being removed 
from NGC 205 by the action of tidal forces, they would 
contribute, at least in projection, to stars with disk-like 
kinematics; we return to this point in §6 below. 

5.5. Field F8 

Figure 21 shows the data for field F8, a field in the pos- 
sible tidal trail of NGC 205. This structure was discussed 
previously by McConnachie et al. (2004b), where they 
presented the morphology of the feature, together with 
kinematics from our field F9 (their field W91). The new 
data in field F8 (adjacent to field F9) do not show the bi- 
modal distribution found by McConnachie et al. (2004b), 
which suffered from a small sample size. Instead the 
only significant stellar population present has, yet again, 
disk-like kinematics with a dispersion of a v = 38 kms -1 
centered at vi ag = — 43kms~ 1 (v'i ag — — 37kms , 



5.6. Field F12 

Figure 22 shows the location of another dense struc- 
ture, this time located towards the North-eastern end of 
the major axis. A peculiarity of this feature (the "North- 
ern Spur" ) is its curious wedge-shaped morphology (seen 
in the upper right-hand diagram) which displays a very 
sharp boundary. The velocity dispersion of the structure 
<7 V = 24kms~ 1 (ay = 23 kms -1 ) is again similar to that 
in other fields. 

5.7. Field F16 

Field F16 is situated in a particularly interesting lo- 
cation, on the North-eastern end of our survey region 
(see Figure 23), where our INT survey detected a very 
extended very low surface brightness structure. It lies 
at the greatest projected radius (46kpc) of all the non- 
stream fields we observed. Zucker et al. (2004) argue 
that this object (which they call Andromeda NE) may 
be one of the most low mass and low surface bright- 
ness galaxies found to date or possibly torn-off mate- 
rial from the disk of M31. The radial velocity distribu- 
tion in the bottom left-hand panel of Figure 23 shows a 
well-defined narrow peak of dispersion a v = 27 kms -1 
centered at vi ag — — 80kms _1 [v'i ag — — 65kms _1 , 
<r v i = 28kms~ 1 ). However, this field is the only one in 
our survey in which there could be some concern about 
Galactic contamination affecting significantly the fitted 
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Fig. 19.— As Figure 16, but for field F6. The Gaussian model 
has viag = 2kms _1 and <r v = 36kms _1 (v'i ag = 4kms _1 , a v i = 
37kms _1 ). The white "L" -shaped hole in the finding chart is due 
to a deliberate shifting of the camera position to centre one of the 
CCDs of the mosaic on NGC 205; part of the inner survey region 
was thereby not covered. The hole in the centre of NGC 205 is a 
consequence of crowding. 



disk-like peak. In Figure 24 we show the Heliocentric 
velocity distribution of stars in this field, with the cut 
at v = — 100 kms" 1 , imposed on all fields to minimize 
the Galactic contaminants. The predictions of the Be- 
sangon Galaxy model are superposed (green histogram) 
normalised to the 16.7' x 5' field of view of DEIMOS, 
but corrected for the ratio of available to observed tar- 
gets in our CMD selection window. Thus according to 
this model, below v = — 100 km s -1 the Galactic contri- 
bution will be minor, and can be neglected. The bottom 
right-hand panel of Figure 23 shows the CMD location of 
the stars within the peak; evidently these are relatively 
metal-rich, and very similar to the other disk-like stars 
found in the fields examined previously. 

The velocity histogram in the bottom left-hand panel 
of Figure 23 shows a spike at velocities — 90kms -1 < 
v iag < — 80kms . The stars partaking in this spike 
are uniformly distributed over the field, and they have a 
similar color spread to the other stars in the sample (the 
corresponding points are circled in the CMD on the bot- 
tom right-hand panel). This lack of distinct properties, 
together with the low statistics indicate that this spike is 
most probably an insignificant (~ 2a) random deviation. 

5.8. Fields F4, F5, F7, F9, F10, Fll, F13, F14 and 
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Fig. 20. — The central panel shows the positions of the targets 
in field F6 with respect to a sketch of NGC 205. The diagonal 
line segments show the minor and major axes of M31. The top 
and bottom panels show the kinematic data as a function of the 
standard coordinate £. The dashed line in the top panel marks the 
systemic velocity of M31, while the dotted line shows the systemic 
velocity of NGC 205. The bottom panel displays the corresponding 
distribution of disk lag velocity, with the dashed line marking the 
locus of circular orbits. 



The remaining outer non-stream fields: F4, F5, F7, 
F9, F10, Fll, F13, F14 and F15, do not individually 
possess a narrow and strong velocity peak near the ex- 
pected velocity of a circularly rotating population (see 
Figure 13). However, by summing the fields together 
we can enhance the signal to detect a global population 
which lags the disk by a similar velocity. The nine fields 
that were added together cover a very large area over our 
survey. Some fields, in particular the halo field F7, are 
almost certainly devoid of any disk stars, so the "disk 
velocity lag" is not a meaningful quantity, and will con- 
tribute only noise to the summed distribution. Never- 
theless, Figure 25 shows that a narrow peak is again de- 
tected, albeit with a large population of "contaminants" , 
this time stemming mostly from the halo population of 
M31 itself. Due to the kinematic substructure in the 
halo population (which will be presented in a forthcom- 
ing paper, Chapman et al., in prep.) we cannot easily 
remove the halo "contamination" to enhance the signal 
of the disk-like peak. Nevertheless, to help guide the eye, 
we have used a maximum likelihood algorithm to fit a 
Gaussian function to the velocity data windowed between 
— 80kms _1 < vi a g < 50kms _1 , which we find to be cen- 
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Fig. 21.— As Figure 16, but for field F8. The Gaussian model 
has vi ag = — 43kms — 1 and tr v = 38kms _1 (v'i ag = — 37kms — 1 , 
<t v i = 35kms _1 ). 



Fig. 22.— As Figure 16, but for field F12. The Gaussian model 
has vi ag = — 78kms _1 and a v = 24kms — 1 (v'i ag = — 60kms — 1 , 
a v i = 23kms~ 1 ). 



tered on v — — 55kms _1 with dispersion a v — 33kms _1 
(v'iag — — 41kms~ 1 , a v > = 31kms _1 ). Although the re- 
liability of this fit is clearly compromised by the strong 
halo contamination, the similarity of this peak with those 
presented in Figures 16 to 23 is very striking. This con- 
firms that a disk-like population is present in almost all 
the surveyed fields. 

5.9. Spectroscopic metallicities 

The spectra also allow a measurement of the metal- 
licity of the targeted stars from the equivalent widths 
(EW) of the Ca II triplet absorption lines. While the 
noise in individual spectra is typically too high to yield 
a useful measurement of the Ca II equivalent widths, we 
proceeded to estimate the average metallicity by stack- 
ing the RGB star spectra of the targets within 2a of 
the fitted peaks of the fields displayed in Figures 16 
to 23 (the corresponding stars are also highlighted in 
the CMD panels). The spectra are shifted to zero ve- 
locity before stacking. We follow as closely as possi- 
ble the method of Rutledge, Hesser & Stetson (1997), 
fitting Moffat functions to the Ca II lines. The aver- 
age spectrum yields a measurement of the Ca II triplet 
equivalent width used to estimate the metallicity as 
[Fe/H] = -2.66 + 0.42[ECa - 0M(V HB - V ave )}, with 
ECa = 0.5EW XS49S + 1.0EW XS542 + 0.6EW XS662 , V HB 
being a surface gravity correction relative to the V mag- 
nitude of the horizontal branch, and V ave the average 
(luminosity weighted) magnitude of the stars in the field. 



We adopt the value of Vhb — 25.17 for M31, measured 
by Holland (1996). The resulting combined spectra are 
displayed in Figure 26, along with the metallicity mea- 
surement (on the Carretta & Gratton 1997 scale). The 
typical uncertainty is ~ 0.2 mags. By chance, the A8498 
line in the summed spectrum in Field 16 falls at the loca- 
tion of a night-sky line, so the two other Ca II were used 
in that instance. For comparison to the inner fields, on 
the top panel of Figure 26 we also display the summed 
spectrum of the disk fields D2-D5, derived from those 
stars within 2a of the fit shown in Figure 10. 

These metallicity measurements are all very similar, 
identical within the uncertainties, and agree reasonably 
well with the position of the targets relative to our fidu- 
cial clusters on the color-magnitude diagram. This yields 
further evidence that the dominant population in the ma- 
jority of the fields we surveyed starting from the disk pop- 
ulations at a radius of 15 kpc out to ~ 50 kpc are probing 
different parts of the same global population, which has 
metallicity of [Fe/H] = -0.9. 

5.10. Scale-length 

A further useful constraint on the nature of this disk- 
like population is its density profile. We proceed by mea- 
suring the number n s of stars within 2a of the narrow 
peaks shown in Figure 10 and Figures 16 to 23, and count 
the number of observed targets n in the CMD selec- 
tion window, and also the number of available targets nt 
within the selection window that fall into the DEIMOS 
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Fig. 23.— As Figure 16, but for field F16. The Gaussian model 
has vi ag = — 80kms _1 and a v = 27kms _1 (v'i ag = — 65kms — 1 , 
<j v i = 28kms" 1 ). 
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Fig. 24. — The heliocentric velocity distribution of field F16, 
situated in the low-surface brightness overdensity to the North-east 
of M31. The green histogram shows the prediction of the Besangon 
model in an area of the DEIMOS field of view, but corrected for the 
fact that only 92 stars in the CMD selection window were observed 
(out of a possible 282). The Besangon model therefore predicts that 
the Galactic contamination is minimal at heliocentric velocities v < 
— 100 km s -1 , which strengthens the results presented in Figure 23. 
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Fig. 25. — The distribution of disk lag velocities for the (non- 
stream) external fields without a prominent narrow disk-like ve- 
locity component; these are: F4, F5, F7, F9, F10, Fll, F13, 
F14 and F15. This combining brings out of the noise, yet again, 
a narrow velocity peak close to expected velocity of a circularly 
rotating population. For illustrative purposes only, we have fit- 
ted the velocity interval — 80kms _1 < Vi ag < 50kms _1 with a 
Gaussian of mean v = — 55kms _1 and dispersion a v = 33kms — 1 
(v'iag = —41 km s — 1 , a v i = 31 km s — 1 ). 
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Fig. 26. — Combined spectra for all the stars within 2a of the 
Gaussian fit. The three Ca II triplet lines at 8498.02 A, 8542.09 A, 

o 

and 8662.14 A are the most obvious features of these spectra. The 
field name is indicated on the left, and derived metallicities on the 
right of each panel. 
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Fig. 27. — The density profile of the extended disk population 
is displayed in the top panel, isolated by counting kincmatically- 
confirmcd stars and scaling by the ratio of possible targets to the 
total number of stars measured in each field. A disk-plane inclina- 
tion of i = 77° is assumed to calculate both the de-projected radius 
R and the de-projected surface density. The dashed line shows an 
exponential fit to the data in the range R < 40 kpc, which implies 
a scale-length of = 6.8 ± 0.4 kpc. The dotted line reproduces 
the fit to the outer major axis profile of Figure 2 (fiR = 5.1 ± 0.1). 
Field F6 is not included in this plot since it lies on the dwarf galaxy 
NGC 205, whose stars enhance enormously the number of possi- 
ble targets in the DEIMOS field. The lower panel shows the same 
data for i = 64.7°. The fitted exponential function has scale-length 
hjj = 6.7±0.3kpc. The dotted line again reproduces the fit of Fig- 
ure 2, corrected for the difference in inclination angle. 



these fields has a slightly different profile to the under- 
lying exponential population revealed in Figure 2, or it 
may also mean that the inclination of the outer disk is 
lower than i = 64.7°. However, it should be recalled that 
the spectroscopic fields preferentially targeted overdensi- 
ties detected in the photometric survey, which need not 
follow closely the medianed major axis profile. Another 
important property to bear in mind when interpreting 
Figure 27 is that the velocity distributions in all fields 
showed narrow peaks at a similar velocity lag. There 
was no indication of bimodal distributions with a com- 
ponent due to an underlying smooth "normal" disk plus 
an additional component due to the sub-structure (with 
the exception of field F2 which we discussed in §5.2). If 
there were an underlying normal disk population it has 
identical kinematics to the substructure. 

5.11. Total luminosity and angular momentum 

In terms of its mass fraction, the outer disk is a very 
significant part of the galaxy. As we have shown in Fig- 
ure 2, the exponential decline can be followed out to 
40 kpc or 7.8 scale-lengths. Approximately 10% of the 
mass of the M31 disk therefore lies beyond the typical 4 
scale- length cutoff (Pohlen, Dettmar & Liitticke 2000). 
This region contains an even larger fraction of the an- 
gular momentum of the disk: assuming a flat rotation 
curve results in ~ 30% of the angular momentum of the 
disk residing beyond 4 scale-lengths. 

5.12. HI content 

Recently, Thilker et al. (2004a) have presented an anal- 
ysis of Green Bank Telescope A21 cm observations of 
M31 and its environment. Their very high sensitivity 
of 2.5 x 10 17 cm~ 2 within a beam-size of 2 kpc allowed 
them to detect many low-mass H I clouds as well as an 
extended filamentary component. A careful inspection of 
their Figure 1, in which the H I data-cube is presented, 
reveals a spatial overlap between the H I emission and 
our extended disk fields. Furthermore, the velocities of 
the H I gas coincide with the peaks of the extended disk 
velocity distributions found in the present study, though 
a detailed comparison is not currently possible as the ve- 
locity intervals in the map presented by Thilker et al. 
(2004a) are large (34kms~ 1 ). 



field of view. The local density in a field is estimated as 
ntn s /n Q — rib, where n& is the background density com- 
puted from fields at the edge of the survey. The top panel 
of Figure 27 shows the resulting RGB star density pro- 
file, under the assumption of the standard Walterbos & 
Kennicutt (1988) inclination angle of i — 11° . The pro- 
file shows an exponential decline to R ~ 40 kpc, at which 
point there appears to be a break in the relation, with 
possibly a flat extension from ~ 40 kpc to ~ 70 kpc. Fit- 
ting the data at R < 40 kpc, we measure an exponential 
scale-length of hn — 6.8 ± 0.4 kpc. This value is some- 
what higher than the 5.1 ± 0.1 kpc scale-length (dotted 
line) that we fitted earlier to the major axis star-counts. 
If we instead assume the lower inclination angle calcu- 
lated in §2 of i = 64.7° (lower panel), the scale-length 
drops marginally to ha = 6.7 ± 0.3 kpc. 

These longer scale-lengths may mean that the popula- 
tion with disk-like kinematics that we have uncovered in 



5.13. Variations in stellar populations 

The color-magnitude diagrams presented in Figures 
16-23 show a very similar broad upper RGB, indicat- 
ing that a similar distribution of stars populates these 
fields. However, recent work by our group (Ferguson et 
al. 2005) using very deep photometry from the Advanced 
Camera for Surveys (ACS) instrument on board the Hub- 
ble Space Telescope reveals significant differences in the 
morphology of the horizontal branch over the inner part 
of the M31 "halo". The surveyed fields include small ar- 
eas of the fields Fl, F5, F9, F10, Fll and F12 from the 
present study. The differences in the horizontal-branch 
populations testify to the existence of inhomogenities in 
the old stellar populations in these fields. 

6. DISCUSSION 

Our radial velocity data reveal the presence of an un- 
expected population, present in almost all the surveyed 
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fields from R = 15kpc out to ~ 70kpc, that lags the 
H I by at most 80kms _1 , and that has a dispersion of 
~ 20kms _1 to ~ 40kms~ 1 . This velocity dispersion is 
large enough to rule out the possibility that these struc- 
tures are bound dwarf galaxies. It should be kept in mind 
that the circularly-rotating disk model that we presented 
in Figure 7 is very simplistic, and one should expect sig- 
nificant deviations from it, especially on fields at large 
y (minor axis distance), where small changes in the in- 
clination of the disk, due for instance to warping, will 
cause large errors in the disk velocity lag. So the good 
agreement in the mean velocity of the disk-like peak in 
each field is all the more surprising. The direction of the 
effect is such that if the disk is warped so that it is seen 
less edge-on, the model of Figure 7 will overestimate the 
rotational velocity, and a population will appear to run 
ahead of circular motion. 

The errors in the model, however, will not affect sig- 
nificantly the fitted <j v in each field. The variations 
in a v in the outer fields are large, ranging between 
<j v = 15kms _1 in field F3 and a v — 38kms _1 in field 
F8. Field F2 (with a v = SSkms -1 ) mirrors field F3 
about the major axis, so we can see that the differences 
in velocity dispersion reflect real changes between fields, 
which arc not just due to different projections of the R 
and (j) components of the velocity ellipsoid onto the line 
of sight. 

For the metallicity measurements we are unfortunately 
forced to co-add the spectra in the velocity peaks to ob- 
tain combined spectra of useful signal to noise. This 
necessarily means that we lose all information on the 
metallicity dispersion of the populations. Nevertheless, 
the mean metallicity can be measured reliably, and the 
derived values in all fields are consistent with the data 
from each field being drawn from a global population of 
identical metallicity. 

All the DEIMOS fields we present in this paper, except 
for fields F5 and F7 were positioned on regions of high 
target density in M31, which is seen as the dark black 
area in a visual inspection of Figure 1. This corresponds 
to a surface number density of stars of approximately 10 4 
degree -2 in the RGB selection window defined in §3. The 
panoramic view of the INT survey provides an invaluable 
complement to the spectroscopic results. Given the very 
similar velocity dispersions, disk velocity lags and metal- 
licities of the stars we have identified, it is clear from 
an inspection of Figure 1 that these populations can be 
united in a single morphological entity: an extended disk- 
like structure which has rather remarkable well-defined 
fragmented sub-structures at its outer edge. 

6.1. Interpretation of velocity lag 

We have found a population that appears to lag the 
velocity of circular orbits in many fields, but what does 
this lag represent? A simple answer may be that, in re- 
ality, the outer disk is less inclined than we modeled. 
However, some fraction of the velocity lag must be in- 
trinsic. In an equilibrium disk, a stellar population with 
non-zero velocity dispersion lags the circular speed; this 
is known as the asymmetric drift (see, e.g., Binney & 
Tremaine 1987, their §4.2). We use Eqn. 4-35 of Binney 
& Tremaine (1987) to estimate the expected asymmetric 
drift, and adopt their simplifying assumptions taken to 
calculate the Solar Neighborhood asymmetric drift re- 
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Fig. 28. — The asymmetric drift of a disk population in equi- 
librium, for three different values of the ^-component of velocity 
dispersion: 20kms~ 1 (dotted line), 30kms _1 (continuous line) 
and 40kms _1 (dashed line). A scale-length of 5.1 ± 0.1 kpc, as 
measured in §2, was adopted. 

lation. Figure 28 shows the expected asymmetric drift 
as a function of radius for disk populations of disper- 
sion n = 20kms~ 1 (dotted line), — 30kms -1 (con- 
tinuous line) and = 40kms _1 (dashed line). These 
estimates show that a significant velocity lag should be 
expected. 

However, there are many uncertainties that hinder a 
direct comparison of these asymmetric drift estimates to 
the data. First of all, the overdensities targeted in the 
outer fields are evidently not equilibrium structures, so 
the asymmetric drift equation need not hold precisely. 
Apart from the uncertainty in the inclination of the 
galaxy, there are further uncertainties in the projection 
of the velocity dispersion ellipsoid, as we discussed in 
§4.1. A realistic treatment of the asymmetric drift will 
require data that sample better the extended disk region. 

6.2. Time-scale for dissolution 

With typical sizes of ^ 5 kpc and velocity dispersions of 
<t„ ~ 30kms -1 , the virial mass associated with the sub- 
structures we discuss in §5 would have to be of the order 
of ~ 10 9 M Q . However, the total luminosity of these 
sub-structures is very low, although this is a difficult 
quantity to measure given the irregular morphology and 
the difficulty of identifying their spatial limits. For the 
case of field F16, where the structure is well-separated 
from M31, the total luminosity is a mere 3 x 10 6 L Q . 
Thus if these structures are bound, they require a ridicu- 
lous mass to light ratio of several hundred in Solar units. 
They are instead most probably unbound structures dis- 
solving in the global potential of the galaxy. 

It is informative to estimate the lifetime of the struc- 
tures if they are indeed unbound. In this case both the 
winding due to differential rotation and the velocity dis- 
persion of the structure set a dissolution time-scale. In 
a potential with a flat rotation curve of velocity v c , an 
unbound structure at galactocentric radius R becomes 
smeared out by an extent equal to its initial size in a 
time R/v c , i.e. in 170 Myr at R = 40 kpc for the case 
of M31 (taking v c = 240 km s -1 ). Alternatively, if we 
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consider just the effect of the expansion of the stellar 
population due to its velocity dispersion unconstrained 
by self-gravity, with a velocity dispersion of ~ 30kms _1 , 
a structure of size 5kpc will double in size in 170 Myr, 
fortuitously the same time estimated before. These sim- 
ple arguments give an order-of-magnitude estimate of the 
the survival time of the structures, and strongly suggest 
that they cannot be long-lived. 

6.3. Age of population 

Brown et al. (2003) present extremely deep photome- 
try with the ACS camera for a field on the minor axis 
of M31. The location of this field, indicated in Figure 4 
with a small square, lies at a projected radius of only 
y = — 11.4kpc. The corresponding de-projected radius 
is R = 53 kpc, beyond the break in the surface density 
profile of Figure 27. However, a tiny 5° warp in the ex- 
tended disk would bring the field into the R < 40 kpc 
region. In their data, Brown et al. find evidence for a 
30% component of intermediate age stars (6-8 Gyr) with 
mctallicity [Fe/H] > —0.5, though they claim a similar- 
quality fit for a population of mctallicity [Fc/H] > —1.0. 
Given that this is the youngest component in that field, 
and has comparable mctallicity, it is possible that the 
intermediate-age population whose age they constrain is 
the extended disk. However, a definitive measure of the 
age of the population will require resolving the main- 
sequence turnoff in fields where we are sure the the ex- 
tended disk population is present. 

6.4. Formation scenarios 

How did the extended disk form? Several possibilities 
present themselves. 

• Formation by the accretion of many small 
sub-galactic structures. This scenario would 
explain naturally the lumpy sub-structure seen at 
the edge of the extended disk, as well as the dif- 
ferences in the velocity dispersion and the varia- 
tions in the old stellar populations between fields. 
However, the chemical homogeneity of the material 
argues against this possibility, since the different 
accreted galaxies would likely have a range of dif- 
ferent mctallicities, depending on their mass and 
star-formation history. Admittedly, the mctallic- 
ity measurements we present above give only the 
mean of the population, and in the central regions 
of M31, considerable mixing must have taken place. 
Nevertheless, it is surprising to find the same mean 
abundance even in the outermost sub-structures. A 
further problem with this scenario is the disk-like 
kinematics. A substantial mass is required for dy- 
namical friction to affect the orbit of a galaxy satel- 
lite (> 1O 8 M , Binney & Tremaine 1987) which 
has to become circularized before being disrupted 
in order to give rise to a final structure with disk- 
like kinematics. Another problem is the disk-like 
distribution of debris. Though accretion in the 
disk plane does appear to be favored (Pcnarrubia, 
Kroupa & Boily 2002), the accretion of a large 
number of halo objects is expected to give rise to 
a distribution closer to that of a spheroidal stellar 
halo than a disk. Even the thick disk component in 
the Abadi et al. (2003) simulation, in which ~ 50% 



is formed by accretion, is a very thick component of 
axis ratios ~ 1 : 3 and scale-length of 4.5 kpc (our 
measurements from their Figures 1 and 3). Form- 
ing the extended disk observed in M31 in this way 
will certainly be a challenge for modelers. 

• Formation from stars accreted in a single 
large merger. Here we envisage the assimila- 
tion of a whole M33-sized galaxy in a single merg- 
ing event. Accreting 10% of the total mass of a 
galaxy counts as a significant merger, affecting the 
morphology of the larger galaxy. However, given 
that the bulk of the interaction took place long 
ago, possibly 6-8 Gyr in the past (corresponding 
to z ~ 1), the progenitor of M31 would have been 
considerably less massive at that time (Hammer 
et al. 2005), and the resulting interaction would 
likely be termed a "major merger" after which An- 
dromeda would have been globally reshaped. Thin 
disks are easily heated up by mergers (Toth & Os- 
triker 1992; Velazquez & White 1999), so this sce- 
nario makes a testable prediction that stars in the 
inner disk of M31 that formed before the merger 
should have a much larger scale-height than the 
younger thin disk. The advantage of this proposi- 
tion is that it accounts naturally for the uniform 
mean abundances measured in the different fields. 
It is questionable that such a merger could have 
given rise to a population with disk-like morphol- 
ogy and kinematics, since typical disk galaxy col- 
lisions result in elliptical end-products (e.g., Hcrn- 
quist 1992). However, new simulations by Springel 
& Hernquist (2005) cast doubt on the inevitability 
of forming elliptical galaxies from such collisions, 
showing that disk galaxies may form even in equal 
mass mergers if the progenitors are gas dominated. 
Further numerical experiments are required to es- 
tablish whether a disk as extended as the one seen 
around M31 can be formed in this way. 

• Formation induced by a single large merger. 

This scenario would be similar to the previous 
one, but the stellar population of the accreted 
galaxy would end up populating the bulge or other 
spheroidal component of M31. This is perhaps con- 
sistent with the presence of a double nucleus in that 
galaxy (Lauer et al. 1993). However the interaction 
itself could have triggered star-formation on a mas- 
sive scale in the pre-existing gaseous disk of M31, 
with gas accreted from the satellite contributing 
to the enrichment. This possibility could explain 
the disk-like structure and kinematics, and has the 
testable prediction that its constituent stars should 
be no older than the interaction itself. The ho- 
mogeneity in mean abundance over the structure 
could be explained by star-formation being limited 
very rapidly before chemical enrichment could take 
place. Such a rapid formation of the disk in this 
way is similar to what is deduced for ~ lmJy sub- 
mm sources at z ~ 2 with star-formation rates of 
~ lOOM /yr (Chapman et al. 2003; Greve et al. 
2005). Nevertheless, this scenario is also problem- 
atic, as it is hard to reconcile with the variations 
in the velocity dispersion between fields. 
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A concern with the last two scenarios is that the sub- 
structures that we see on the outskirts of the extended 
disk now are short-lived features that will soon disappear 
(<< lGyr). This means that we would be observing 
the system at a special time at the very tail-end of the 
merging process. A further problem is the constraint on 
the size of a past major merger set by the existence of a 
population of thin disk globular clusters that has recently 
been discovered in M31 by Morrison et al. (2004). Those 
authors infer that no more than 10% of the mass of the 
disk could have been accreted in a single merger since 
the formation of the clusters. The time of formation of 
the clusters is currently poorly constrained, though they 
find no reason why the clusters should be significantly 
younger than those of the Milky Way. Measuring the 
ages of these objects will help greatly in distinguishing 
between the above scenarios. 

6.5. Possible connection to NGC 205? 

As we noted in §5.4, if stars become unbound from 
NGC 205, we would observe them in field F6 with kine- 
matic properties very close to that of disk material. Like- 
wise, if McConnachic ct al. (2004b) were correct in iden- 
tifying the arc- like structure to the North of M31 as a 
tidal stream from NGC 205, this tidal debris, which is 
located in fields F8 and F9, would also have disk-like 
kinematics. 

This raises the possibility that NGC 205 is connected 
to the extended disk structure, cither as the direct sup- 
plier of stars, or as the remnant of the ancient structure 
which created the extended disk. It is interesting to note 
in this context that the mean mctallicity of NGC 205 has 
been measured to be [Fe/H] ~ —0.9 (Mould, Kristian & 
Da Costa 1984), identical to the mean spectroscopic mea- 
surements from the fields presented in this survey. 

A counter-argument however, is that recent estimates 
for the distance to NGC 205 indicate that it is behind 
M31, being at a distance of 824 ± 27kpc (McConnachie 
et al. 2005), compared to 785±25-40.6sin(i) kpc for the 
M31 disk at the location of NGC 205. The difference in 
distance is therefore 79 ±37 kpc, little over 2a, although 
we are probably grossly overestimating the uncertainties, 
since a differential distance measurement is expected to 
be much more accurate. 

However, this need not rule out the connection with 
NGC 205. If the dwarf galaxy were indeed ~ 80 kpc 
beyond the foreground disk on the near side of M31, it 
would be <~ 40 kpc behind M31, at a location ~ 17 kpc 
above the plane of the disk on the far side of that galaxy. 

6.6. Comparison to Hurley-Keller et al. (2004) 

Hurley-Keller et al. (2004) have presented an analysis 
of the kinematics of 135 planetary nebulae in the disk and 
inner halo of M31. They find the surprising result that 
a standard model for the thin disk, thick disk and bulge 
cannot reproduce the observed PNe kinematics at major 
axis distances of 15 kpc to 25 kpc. In order not to intro- 
duce the presence of an additional component, they in- 
voke Occam's razor to alter the parameters of the bulge, 
deducing that the bulge must rotate rapidly at large dis- 
tances. We re-assess their data in the light of the discov- 
ery of the additional extended disk population. Figure 29 
shows the positions of the 135 PNe in their survey. Stars 
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Fig. 29. — The location of the PNe from Hurley-Keller et al. 
(2004) with respect to M31; the lines and ellipses from Figures 1 
and 4 are reproduced to show the scale. The small open circles 
show the full sample of 135 PNe, while our chosen sub-sample of 
77 objects with \x\ > 8 kpc are indicated with filled circles. 



close to the centre of M31 are a complex mix of popula- 
tions, so we investigate the sub-sample with |a;| > 8 kpc, 
which are displayed as filled circles. The distribution of 
velocity lags is shown in Figure 30. We fit this distribu- 
tion for vi a g > — 150 km s -1 using a maximum-likelihood 
algorithm, to account for the velocity measurement un- 
certainties quoted by Hurley-Keller et al. (2004) as well 
as the disk model uncertainties (from Figure 8). The 
resulting Gaussian model has vi ag — — 52kms~ 1 and 
a v = 44kms _1 , very similar to the fit in Figure 10. 
Further analysis is necessary to establish conclusively 
whether the existence of the extended disk alleviates the 
need for a rapidly rotating outer bulge as suggested by 
Hurley-Keller et al. (2004). 



6.7. Globular clusters 

A population of massive globular clusters has recently 
been found in the disk of M31 (Fusi Pecci et al. 2005). 
These objects are scattered around the outskirts of the 
inner disk at distances of ~ 20 kpc, and have very 
blue colors, with color indices indicating a young age 
of ~ 2 Gyr. The population is also very significant, ac- 
counting for £ 15% of the globular cluster system of 
that galaxy. This population of young clusters rotates 
rapidly, possessing kinematics similar to that of the disk. 
A further piece of relevant information is provided by 
the recent discovery of the globular cluster Bologna 514 
(Huxor et al. 2003; Galleti et al. 2005), which is situ- 
ated close to the South-western major axis of M31 at a 
projected distance of 55 kpc, almost exactly on the outer 
ellipse marked in Figures 1 and 4. The rotation model of 
Figure 7 gives a disk-lag velocity of this globular cluster 
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Fig. 30. — The distribution of disk lag velocities of the Hurley- 
Keller et al. (2004) PNe with |x| > 8kpc. A maximum-likelihood 
Gaussian fit to the data in the interval — 150 km s -1 < V[ ag < 
100 km s -1 is also shown, for which both the estimated velocity 
uncertainties from Hurley-Keller ct al. (2004) as well as the un- 
certainties in the disk rotation model (from Figure 8) have been 
taken into account. We fit a Gaussian mean velocity lag of vi ag = 
— 52kms~ 1 and dispersion <r v = 44kms~ 1 (v'i ag = — 44kms~ 1 , 
a v i = 46kms~ 1 ). The same low dispersion disk-like population as 
discovered in our DEIMOS fields is also present here. 



of vi a g = -70kms _1 (v'i ag = -53kms _1 ). 

Though it is tempting to speculate on a possible con- 
nection between these globular clusters and the extended 
disk, the metallicities of the clusters appear to be much 
lower than the [Fe/H] = —0.9 we measure for the struc- 
tures in this study. However, as pointed out by Fusi Pecci 
et al. (2005), the metallicities of these globular clusters 
could be higher if they are of younger age than the typical 
halo cluster. Further observations are needed to estab- 
lish the age of these M31 disk clusters, and their relation 
to the extended disk. 

6.8. A connection with the thick disk? 

Several edge-on galaxies have been found to possess 
a flattened component with a much larger scale-height 
than the thin disk (Dalcanton & Bernstein 2002) proba- 
bly indicating the presence of a thick disk. Recently, two 
of the sample of Dalcanton & Bernstein (2002) were ob- 
served spectroscopically by Yoachim & Dalcanton (2005), 
who found the thick disk to be rotating very slowing in 
one case and counter-rotating in the other. This result 
shows that thick disks can have completely different kine- 
matics from the thin disk, although this is not the case 
for the Milky Way. The fact that the extended disk we 
find around M31 has a low velocity dispersion, and fol- 
lows circular orbits closely thus behaving kinematically 
like a thin disk, as well as having a very similar scale- 
length to the thin disk, suggests a close link with that 
component, and argues against association with the thick 
disk of Andromeda. It should be noted, however, that 
it has not yet been shown that Andromeda does indeed 
possess a thick disk. The present dataset may allow us 
to detect the thick disk, but we will return to this issue 
in a future contribution. 



6.9. Extended disks in other galaxies? 

Our group has also been attempting to unravel the 
formation history of the Milky Way (Ibata et al. 2001a, 
2002), with similar goals to our study of M31. Following 
up on the discovery of a peculiar halo structure in the 
SDSS (Newberg et al. 2002), we showed that the Milky 
Way is enveloped in a Ring-like distribution of old stars 
that lie at Galactocentric distances of <~ 20 kpc, well be- 
yond the edge of the Galactic disk (Ibata et al. 2003). A 
subsequent analysis of asymmetries in the 2MASS sur- 
vey revealed a strong enhancement of stars at the edge 
of the Galactic disk, which we interpreted as an accreted 
dwarf galaxy fragment (Martin et al. 2004a, b; Bcllazzini 
et al. 2004, see also Rocha-Pinto et al. 2005); a connec- 
tion with the ring-like structure is possible, though not 
proven. However, this "Galactic Ring" or "Galactic An- 
ticenter Stellar Structure" (Crane et al. 2003) is not the 
only stellar population with disk-like kinematics that lies 
beyond the rim of the Milky Way thin-disk: recently Ma- 
jewski et al. (2004) and Rocha-Pinto et al. (2004) have 
uncovered a further stellar structure located at a dis- 
tance of ~ 40 kpc and which spans a huge angle on the 
sky. These findings indicate that the outer regions of the 
Galactic plane have a similar messy flattened structure 
to what is observed in M31 (Figure 1). 

A possible young counterpart to the "extended disks" 
in the Milky Way and M31 has been seen in some more 
distant galaxies. GALEX observations of M83 (Thilker 
et al. 2004b) have revealed an extended population of 
UV sources (OB stars) beyond the inner disk (though 
at much smaller radius than the structure we detect in 
M31), while several galaxies are known to exhibit H II 
regions to extremely large radii (Ferguson et al. 1998). 
Thus it appears that there may exist an additional com- 
ponent to spiral galaxies, a putative "extended-disk" . 
The implications for galaxy formation arc profound: the 
inhomogenous spatial distribution and moderately high 
metallicity suggests an origin in multiple accretions, yet 
accreted fragments, which are the building blocks of 
galaxies, do not give rise to these huge extended disks 
in current simulations. 

7. CONCLUSIONS 

We have presented an analysis of a large survey of 
the Andromeda galaxy with the Keck/DEIMOS spectro- 
graph based on targets selected from a panoramic survey 
with the INT Wide Field Camera. A huge rotating struc- 
ture that follows an exponential density law with scale- 
length 5.1 ±0.1 kpc is detected out to a de-projected dis- 
tance of ~ 40 kpc, with sub-structures continuing out 
even further in radius. All of the photometric sub- 
structures that were targeted for spectroscopic follow-up 
participate in the same disk- like movement. The velocity 
dispersion is low, typically 30kms _1 . These properties 
indicate that the population is a huge extended stellar 
disk. The color dispersion on the RGB of stars that 
belong to this structure is broad, implying that a wide 
metallicity range is present. Spectroscopic measurements 
of the mean metallicity give [Fe/H] = -0.9 ± 0.2. No sig- 
nificant radial trend in cither kinematics or metallicity 
is found, although deep ACS photometry in several of 
the DEIMOS fields reveals variations in the horizontal 
branch morphology over the galaxy, indicating that the 
stellar population of the most ancient stars differ between 
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fields. We estimate that ~ 10% of the luminous mass of 
the Andromeda disk resides in this extended disk compo- 
nent, beyond the 4 scale- length cutoff observed in most 
galaxies. However, this region may account for ~ 30% of 
the total angular momentum of the galaxy. 

The spatial sub-structures that form the outer edge 
of this system are almost certainly no longer gravita- 
tionally bound and will disappear as spatially coherent 
entities within ~ 200 Myr. However, the existence of 
these sub-structures indicates a formation mechanism by 
mergers, and implies that the extended disk is still form- 
ing, although we seem to be observing it at the very tail 
end of its formation period. It is unclear at this point 
whether the extended disk was formed by multiple accre- 
tion events or by a single large event, though the latter 
possibility appears more plausible given the present data. 

Many questions remain to be answered. Is the ex- 
tended disk a continuation of the thin disk, or is it a 
separate component? Was the merger, or were the merg- 
ers, that formed the extended disk responsible for the for- 
mation of the i? 1 / 4 -law halo of Andromeda? Why is the 
extended disk closely aligned with the thin disk? Why is 
it not extremely warped? 
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The surface brightness of this structure, 27- 
31 mag/arcsec 2 in its outer regions, places it beyond 
what would be classified as "disk" in more distant sys- 
tems. Though such faint surface brightness levels are 
extremely challenging to observe, by resolving individ- 
ual stars it will be possible to ascertain whether the ex- 
tended disk is a particularity of M31 or whether similar 
components are present in other galaxies. Further work 
is also needed in the Milky Way to study the distant 
stellar populations hidden by foreground populations in 
the Galactic plane, and to understand the nature of the 
structures recently unveiled at the edge of the Galactic 
disk and beyond (Newberg et al. 2002; Ibata et al. 2003; 
Martin et al. 2004a,b; Majewski et al. 2004; Rocha-Pinto 
et al. 2004, 2005), and their relation to the vast extended 
disk in the Andromeda galaxy. 
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